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The low frequency quantum transport properties of a three-probe mesoscopic 
conductor are studied using Biittiker's AC transport formalism. The static trans- 
mission coefficients and emittance matrix of the system were computed by explic- 
itly evaluating the various partial density of states (PDOS). We have investigated 
the finite size effect of the scattering volume on the global PDOS. By increasing 
the scattering volume we observed a gradual improvement in the agreement of 
the total DOS as computed externally or locally. Our numerical data permits a 
particular fitting form of the finite size effect. 
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1 Introduction 



In recent years extensive investigations of ballistic and mesoscopic quantum con- 
ductors have been carried out[|IJ. Experimentally the advances in nanotechnology 
have enabled the possibility of fabricating submicron structures with linear size 
of 1000 Angstrom or less. Due to quantum size effect, the transport properties 
of these small systems can be very different from their classical counterpart and 
many interesting phenomena have been discovered On the theoretical side, a 
main tool for understanding static ballistic transport is based on the scattering 
approach of Landauer-Buttiker formalism ||. [|. 

It has been realized that the usual static scattering approach can not be 
directly applied to dynamic transport problems where the external potential has a 
time dependent oscillating component. As shown by Biittiker and his co-workers, 
a direct application of the original approach of Landauer-Buttiker formalism can 
not yield electric current and charge conservation. To preserve this conservation, 
it is necessary to consider the implication of the long range Coulomb interaction. 
As a result the AC transport theory for coherent quantum conductors are more 
complicated. At present there are several approaches to deal with the problem 
of computing AC conductance. In a strongly correlated electronic system, the 
Anderson impurity model is often used. To treat AC transport one employs the 
linear response theory in conjunction with the Keldysh Green's function which is 
often applied to deal with non-equilibrium problems^. One can also use Kubo's 
linear response theory by assuming that the electric field inside the sample is 
known a priori. However this is a very strong requirement []5j . 

Along another line of development Biittiker and his co-workers have advanced 
a current conserving formalism ||. The key idea in this theory is to consider the 
self-consistent internal potential so that the current and the charge is conserved. 
In a series of articles, Biittiker and co-workers investigated several low frequency 
quantum transport problems || |7|, |[ §, [it], |TIJ. To first order in frequency u, the 
response of an arbitrary scattering problem to quasi-static perturbations in the 
scattering potential is naturally expressed in terms of a set of local partial density 
of states (PDOS) each associated with one element of the scattering matrix. This 
AC transport formalism has also been extended to second order in frequency uj 



in the quantum hall regime 11 



The application of Biittiker's AC transport theory is easier in ID, such as 
a ID quantum well, a 5-function potential, and a perfect quantum wire, where 
the scattering matrix and wave function can be obtained analytically. Much 
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intuition and interesting results have been obtained from ID calculations which 
can often been done analytically. The investigations for 2D conductors have 
recently begun [12, [13|, [14], FjJ and the dynamic transport response of a 3D atomic 



wire has been calculated using first principles [[If]]. However, due to technical 
difficulties, in 2D one usually can not obtain analytical expressions for the AC 



admittance except for very special cases [13, 14]. For more general situations 
numerical calculations are needed. 

It is the purpose of this article to further investigate AC admittance of 2D 
coherent conductors in the ballistic regime. In particular we shall focus on numer- 
ically analyzing a three probe conductor as shown in Figure (1). There are several 
motivations for this study. First, similar to that of the DC transport situation^, 
we believe coherent AC admittance of multi-probe (by multi-probe we mean more 
than two probes) conductors should be studied in detail since usually experimen- 
tal measurements are conducted in multi-probe setups. However to the best of our 
knowledge there is as yet a detailed numerical analysis of any 2D multi-probe sys- 
tems. Second, in our investigations of two-probe conductors []l"2"|, |13|, an important 
technical point is the size effect of the scattering volume. It was found |18|, [12], [13|] 
that the total PDOS as computed from external global PDOS (GPDOS) does 
not equal to that computed from the local PDOS (LPDOS), unless the scattering 
volume is very large. This led to a violation of the current conservation and 
gauge invariance in numerical calculations where the scattering volume is always 
finite. Hence there is a need to explicitly and systematically examine the scat- 
tering volume size effect. Finally, in order to study certain physical effects such 
as the inelastic and dissipative effects using the quantum scattering approach, a 
very useful phenomenological approach is to introduce fictitious links from the 



conductor to external dissipative reservoirs [17|]. In this case one must deal with 
multi-probe situations. 

To these purposes we have computed the low frequency admittance of a three- 
probe junction (see Fig. (1)). We have examined the emittance matrix for both 
the tunneling regime and the transmissive regime in detail. The behavior of the 
emittance matrix is found to be closely related to that of the transmission coeffi- 
cients. We have computed the total PDOS from both the external GPDOS and 
the local PDOS, and largely speaking the two total PDOS as obtained approach 
each other as the scattering volume is increased. We found that there exists a 
"critical region" in energy near the second propagating subband threshold, such 
that within this region the charge conservation is not strictly obeyed for any fi- 
nite scattering volume. However the larger the scattering volume, the smaller the 
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"critical region" is. 

This paper is organized as the following. In the next section we set out the 
theoretical and numerical procedures of computing the transmission functions. 
In section III we present and analyze the numerical results. Both the dynamic 
and static transport properties and their relationship will be discussed. Finally, 
a summary is given in section IV. 



2 Theoretical and numerical analysis 



The current conserving dynamic transport formalism proposed by Biittiker and 
co-workers is amply reviewed in several articles^] and we refer interested readers 
to them. In this section we shall only outline our theoretical and numerical 
procedures to compute the necessary quantities such as the various partial density 
of states for the 3-probe system. 

It has been shown by Biittiker, Thomas, Pretre, Gasparian and Christen 0, 
|6|, [7], [18| , to linear order in frequency u the admittance is given by the following 
formula, 

<4M = 9 e a ^ = 0) - tue 2 E a p (1) 

where the emittance matrix E a p is calculated from the various partial density of 
states: 



E 



dN, 
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dE 
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(2) 



The subscripts a/3 indicates scattering from a lead labeled by f3 to that labeled 
by a. The first term in the emittance matrix gives the AC response of the sys- 
tem to the external potential change, while the second term is from the internal 
potential change induced by the external perturbations. The external contribu- 
tion is determined by the global partial density of states |J: for a large scattering 
volume the global PDOS can be expressed in terms of the energy derivative of 
the scattering matrix elements! 
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On the other hand the internal contribution D a p is related to the local PDOS 
and within the Thomas-Fermi linear screening model is given by 
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where the local PDOS (called injectivity) is calculable from the scattering wave- 
function, 

where J is the incident flux and ^ Q (r) is the scattering wavefunction for electrons 
coming from the probe a. In the absence of a magnetic field, the emissivity 



dn(/3,r)/dE equals to the injectivity ||TT | . Finally, dn(r)/dE = Y. a dn(a 1 r)/dE is 



the total local density of states. It is straightforward to prove that the current is 
conserved since the admittance matrix g 1 satisfies S Q (^(cu) = 0. This can be 
seen by realizing that J2a dN a p/dE = dNp/dE is the injectance which is identical 

to J2a D a (3. 

To compute the various PDOS, for simplicity we shall focus on the first trans- 
port subband only, thus the incoming electron energy is restricted to within the 
interval (vr/a) 2 < E < (2-7r/a) 2 in units of h 2 /(2ma 2 ) with m the effective mass 
of the electron and a the width of the leads (see Fig. (1)). Multi-subbands can 
also be included without difficulties, such as that of Ref. [PHI . The scattering 
properties of the three-probe system is then characterized by a 3 x 3 scattering 
matrix S(E) = {s Q/ g} with a, (3 — 1,2,3. For example, for an incident electron 
coming from probe 1, it scatters in the scattering volume, and then reflects back 
to probe 1 with a probability amplitude given by |sn|, or transmits to probes 2 
and 3 with probability amplitudes | S21 1 and |s 3 i|, respectively. The transmission 
coefficients can thus be expressed in terms of scattering matrix, i.e. T a p = |s a( a| 2 . 
For the system of figure (1), the scattering matrix has the following symmetry: 
= I s 2 2 1 , |s 2 i| = I S12 1 , I s 3 i I = |s 32 |, |s 13 | = |s 23 | and |si 3 | = |s 3 i|. Therefore, 
there are only four distinct elements out of nine. 

For the three-probe conductor of Figure (1), the quantum scattering problem 
is solved using a mode matching method. The wavefunction in region I can be 
written as 

*/ = Ex»(y)Ke*" + 6»e" <M ) ) (6) 

n 

where Xn{y) is the transverse wave function, k\ = E — (nrr/a) 2 is the transport 
energy, a n is the input parameter, and b n is the reflection amplitude. Similarly 
for region II, we have 

*!/ = E Xn(y)(c n e th " x + d n e~^ x ) . (7) 

n 

For region III, 

*in = J2Xn(x)(e n e lk -y + f n e~ ik ^), (8) 



5 



where c n and e n are transmission amplitudes and d n and /„ are input parameters. 
The wavefunction in region IV is the combination of wavefunctions in regions I, 
II, and III. At the boundaries of the various regions, we match the wavefunctions 
and their derivatives and this leads to the desired transmission coefficients with 
which the scattering wave functions Eqs. are also determined. 

If we choose point O as the origin (see figure 1), the scattering matrix s\p is 
defined as 

sn = h 

sis = e ie * a . (9) 

Using Eqs. (@)-(0) and the solution of the scattering problem, we can explicitly 
compute the low frequency admittance. 



3 Results 

We have investigated the transmission coefficient and the emittance matrix in 
two different transport regimes for various system parameters. The first regime 
is very transmissive and the second is a tunneling regime where tunnel barriers 
are added at the probes. The AC response of these regimes can be quite different 
as a transmissive situation tends to be inductive, while a non-transmissive case 
tends to be capacitive (see below). The low frequency admittance is given by 
Eq. (|l|) in which the DC conductance g^pioj = 0) of our three-probe system is 
determined using transmission coefficients by applying the Biittiker multi-probe 
conductance formula ||. 



3.1 Emittance 

In Figure (2) we show the transmission coefficients and the emittance E a p in the 
transmissive regime as a function of the incoming electron energy. In this case the 
system does not show any resonance behavior and the transmission coefficients 
T a p(E) are quite large for most of the energy range while the reflection coefficient 
R\x is small (Fig. (2a)). It is interesting to find that the shape of emittance 
are similar to that of the corresponding transmission coefficients, as shown in 
Fig. (2b). This is different from cases where quantum resonances dominant 
the transport |L2) (see below) and for that case the AC responses follow the DC 



transmissions only at the resonances. 



6 



There are two different responses to the external time varying potential: ca- 
pacitive and inductive depending on the sign of the emittance matrix element 
En. According to Eq.(|2|), En consists of two terms: dNu/dE the capacitive 
term and Dn the inductive term. For a two probe capacitor there is no DC cur- 
rent so that dNi 2 /dE = 0. As a result Ei 2 is negative. Therefore for a capacitor 
En = —E12 is positive. Extending this notion, one concludes that the system 
responds capacitively if En is positive. For a ballistic conductor with complete 
transmission dNn/dE vanishes and En is negative. In other words, negative En 
gives an inductive response. These different responses are clearly shown in Fig. 
(2). 

The AC transport properties are very different in the tunneling regime. To 
establish such a regime, we have put tunneling barriers inside probes 1 and 2 
at the junctions between the probes and the scattering volume. In particular 
the barrier heights are Vbarrier = 40-Ei, and the width is 0.1 where the width of 
the wire a has been set to one. No barrier is added in probe 3. We have also 
included a potential well with depth V we u = —40Ei in the center of the scattering 
volume with a size of 2.8 x 1.9. The well and barriers establish several transport 
resonances, these are clearly marked by the sharp peaks in the electron dwell time 



defined as 21 



r a = -l |vMr)| 2 d 3 r , (10) 



J 

where Q is the scattering volume. T\ is plotted against energy in Fig. (3) while the 
inset shows T3. The dwell time measures the duration an electron spends in the 
scattering volume. Thus if transport is mediated by resonance states we expect 
much longer dwell times W\ at the resonances. This idea has recently been proved 
by Iannaconne p3f . Fig. (3) shows that three resonance states, with energies Ei = 
13.2, E 2 = 24.1 and E s = 35.6 are established. The quantum resonances also 
leads to sharp peaks in the transmission coefficient T 2 i and reflection coefficient 
Rn, as shown by the solid lines of Fig. (4a,4b). At these resonances both the 
GPDOS and LPDOS take maximum values, leading to the sharp jumps in the 
emittance En and E 2 \ as shown by the data points in Fig. (4a,4b). The variations 
of En and E 2 \ as functions of energy E are very closely correlated with those of 
Rn(E) and T 2 i{E) near the resonances, as shown by Fig. (4). Since there was no 
tunneling barrier in probe 3, the resonance transmission to that probe is not as 
sharp, and the transport behavior shows a mixture of tunneling and transmissive 
nature, as shown in Fig. (4c). 

In the tunneling regime the AC response changes sharply from inductive at 
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one side of the resonance energy to capacitive on the other side of the resonance or 
vise versa, in distinctive difference as compared to the transmissive case discussed 
above. Let's examine En near resonance E 3 . As the energy approaches to E 3 , the 
system first responds inductively and is followed by a strong capacitive response. 
This behavior is clearly related to the fact that the resonance is characterized by 
a complete reflection indicated by the large peak in the reflection coefficient (see 
Fig. (4a)). This behavior has been seen previously in 2D quantum wires El. On 
the other hand, for ID resonance tunneling, a Breit- Wigner type transmission 
resonance gives rise to the similar AC response behavior discussed here. When 
the incident energy is near the resonance Ei, the AC response is reversed: first 
capacitively and then inductively. Hence the behavior near E\ and E% are very 
different. For energy near E\ the emittance behaves like an odd function but 
near E% it is like an even function. The reason, as we have checked numerically, 
is that the external and the internal responses do (not) reach the maximum at 
the same energy for E near E 3 (Ei). This behavior of E n is also a manifesta- 
tion of the reflection coefficient R\\. As energy sweeps through E%, the strong 
capacitive AC response is due to the complete reflection peak, and the following 
inductive response is because the reflection coefficient Ru « 0. Hence in the AC 
response of a system, near a quantum resonance whether it is voltage following 
current (capacitive) first, or current following voltage (inductive) first, can only 
be determined by detailed analysis and the outcome depends on the peculiarities 
of the system such as the existence of a third probe as we have studied here. 
In Figure (4c) we show the emittance matrix elements E13. Although they have 
much smaller values they do exhibit dips around three resonant energies Ei, E 2 , 
and E3. 

3.2 Finite size effect of the GPDOS 

A very important formal advance of the AC transport theory is the correct char- 
acterization of electric current conservation. In principle this requirement is sat- 
isfied by the AC transport formalism used here|| which demands J2 a Ea/3 = 0. 
Hence we must have J^adN^/dE = J2aD a f3- Since both sides of this equation 
represent the total scattering DOS, thus the current conservation is obtained. In 
practical calculations, the left hand side of this equation is computed externally, 
using the scattering matrix which is calculated at the boundaries of the scattering 
volume. On the other hand the right side of this equation is calculated locally, 
using the scattering wavefunction inside the scattering volume. These two quan- 
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tities becomes precisely equal when the scattering volume is very large[19]. For 
a finite scattering volume, correction terms should be added to the GPDOS, as 
shown in Ref. p|] for ID systems, and in Ref. |13[ for a 2D system. Without 



the corrections, numerical results for a 2D quantum wire showed a small but 
systematic deviation from the precise current conservation [|T2|j . Such a deviation 
actually diverges at the edges of successive propagation subbands as shown in 
Ref. 0. 

Since partial density of states play a vital role in the AC transport formalism 
used in this work, in this section we present a detailed analysis of the finite size 
effect of the scattering volume to the GPDOS. To this purpose we have examined 
a variety of system sizes L for many energies near the onset of the second transport 
subband. As a measure, we define a quantity which is the difference of the total 
DOS as calculated from GPDOS and LPDOS: 

• (id 

Obviously 5p = if the current is precisely conserved. 

Fig. (5a) shows Si as a function of the system size L for three energies very 
close to the second subband edge which is located at £2 = 39.4784. A clear 
crossover to the large volume limit is revealed as 5\ — > when L is increased. It 
is also clear that for energy closer to £2, the crossover is slower (solid line). We 
found that the decay of 5\ is essentially exponential for all energies examined, 
and has a interesting form for large L: 

S 1 ~ e~ fe2(2L+1) (12) 

where (2L + 1) is precisely the scattering volume length from probe I to probe II, 
and k 2 is the momentum corresponding to the second subband energy £ 2 . We have 
plotted — ln[5x/(2L+ l)]/k 2 in Fig. (5b) for several energies. Our numerical data 
supports Eq. ( |T2"D quite well for large L, and for energies closer to £ 2 . It is not 
difficult to understand the form of Eq. flT2"|). Due to the scattering at the junction 
where the three probes meet, complicated mode mixing takes place. While the 
incoming electron is in the first subband, mode mixing generates wavefunctions for 
many higher subbands, including the second subband, which become evanescent 
in the probes. For a scattering volume with a small L, the evanescent mode 
may "leak" out of the volume. However when we calculate the GPDOS from the 
scattering matrix, these "leaked" evanescent modes are not explicitly included, 
leading to a finite 5p. As we increase L, the evanescent modes decays away, and 
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Sff is reduced. In a specific example which can be solved exactly [jl3fl , a similar 
form to Eg. (|T2|) was derived which was needed to correct the GPDOS in order 
to satisfy the precise current conservation. Our numerical study presented here 



reinforces the results of Ref. |T3 



To further investigate the finite size effect to GPDOS, in Fig. (6a,b) we plot 
the total DOS as obtained by GPDOS and LPDOS as functions of energy, for 
three system sizes L. The current conservation condition is satisfied very well for 
most of the first subband energies. When approaching the end of first subband, 
the current conservation condition is violated gradually, i.e. 5p ^ 0. We see 
that for the smallest scattering region L = 0, the agreement of the two total 
DOS is at best reasonable when the incident electron is from probe I and is away 
from the second subband edge (Fig. (6a)), and is quite bad when the electron is 
coming from probe III (Fig. (6b)). The situation improves considerablely when 
we increased the system size. As shown in Fig. (6), for L = 1 and L = 2, the 
agreement of the two total DOS are much better. However there is always a 
divergent behavior near the second subband for all sizes examined if the energy 
is made very close enough to £2- Hence the effect of increasing the size of the 
scattering volume is to decrease the "critical region" where the two total DOS 
disagrees. 



4 Summary 

In conclusion, the low frequency quantum transport properties of a three-probe 
mesoscopic conductor are studied using Buttiker's current conservation formal- 
ism. The static transmission coefficients and emittance matrix of the system 
with different electric potentials are computed. We found that the behavior of 
the emittance matrix is closely related to that of the transmission coefficients. 
We examined the finite size effect of the GPDOS which affects the electric cur- 
rent conservation. In general as the incoming electron energy E approaching the 
threshold of the second subband, the finite-size GPDOS diverges and the current 
conservation is violated. The effect of increasing the size of the scattering volume 
is to decrease the region where the current conservation is violated. 
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Figure Captions 



Figure 1. Schematic plot of the 3-probe quantum wire system. The scattering volume 
is defined by the dotted lines. 

Figure 2. The transmission coefficients and the emittance E a p as functions of the 
incoming electron energy without the tunneling barriers, (a). Solid line: 
reflection coefficient Ru; dotted line: transmission coefficient T 2 i] dashed 
line: T 31 . (b). Emittance E a @. Solid line: E\\\ dotted line: E 21 ; dashed 
line: E$\. 

Figure 3. Electron dwell time t\ as a function of the incoming electron energy in the 
tunneling regime. The three peaks indicate three resonance states in the 
system in this energy range. Inset: T3. 

Figure 4. The transmission coefficients and the emittance E a p as functions of the 
incoming electron energy in the tunneling regime, (a). Solid line: R\\] 
dotted line: En; (b). Solid line: T 2 \\ dotted line: E 2 i] (c). Solid line: T31; 
dotted line: E 31 . 

Figure 5. (a). The difference, Si, of the total PDOS as computed from the GPDOS 
and LPDOS from Eq. ([11]) as a function of the scattering volume linear size 
L. Solid line: at energy E = 39.46699; dotted line: at energy E = 39.45371; 
dashed line: at energy E = 39.44137; (b). The quantity, —ln[5i/{2L+l)]/k 2 
as a function of the linear size L for several incoming electron energies as 
shown. At large L, this quantity approaches unity, confirming the form of 
Eq. (O). 

Figure 6. Comparison of the total PDOS computed from the GPDOS and LPDOS, as 
a function of the incoming electron energy for three different sizes L = 0,1, 
and 2 in the transmissive regime, (a). Electrons coming from probe 1. (b). 
Electrons coming from probe 3. The agreement of the total PDOS is quite 
good up to the "critical region" near the onset of the second subband. 
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TempleGraph 2.4, Origin: /morgoth/users/guo@celeborn - Tue Jan 21 10:29:01 1997 (data file was fig5b.dat) 
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